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Abstract 

We study statistical distributions in a mechanical model for an earthquake 
fault introduced by Burridge and Knopoff [R. Burridge and L. Knopoff, Bull. 
Seismol. Soc. Am. 57, 341 (1967)]. Our investigations on the size (moment), 
time duration and number of blocks involved in an event show that exponen- 
tial distributions are found in a given range of the paramenter space. This 
occurs when the two kinds of springs present in the model have the same, 
or approximately the same, value for the elastic constants. Exponential dis- 
tributions have also been seen recently in an experimental system to model 
earthquake-like dynamics [M. A. Rubio and J. Galeano, Phys. Rev. E 50, 
1000 (1994)]. 
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I. INTRODUCTION 



Systems that present stick-slip dynamics and scaling invariance have attracted consider- 
able attention recently This was triggered, in part, by a seminal paper of Bak, Tang 
and Wiesenfeld, which showed that a class of systems presenting avalanches ou earthquake- 
like dynamics can naturally attain a critical state characterized by power-law distributions 
||T|]. They denoted this phenomenon Self- Organized Criticality. 

Recently, two experimental studies in a continuous elastic system have been performed 
to model the stick-slip dynamics observed in earthquake- like phenomena. In one of them , 
a glass rod is pulled on top of a latex membrane. In this system the slipping events appear 
as detachment waves and a wide distribution of event sizes was observed. However, robust 
scaling behavior was not seen. 

In the other system Q, the elastic medium consists of a transparent gel sheared between 
two coaxial circular cylinders. The inner cylinder is rotated at very low angular speed. Thus, 
this system closely resembles the physics of a spring-block model 0] with periodic boundary 
conditions. Due to the friction between the gel and the cylinders, stick-slip behavior is seen 
in the form of detachment waves. Rubio and Galeano were able to identify four regimes in 
this system, by varying the gel's rigidity, the rotor angular speed and the friction properties 
of the gel with the inner cylinder. The regimes identified were: (a) uniform slipping of the 
gel with respect to the rotor; (b) relaxation events, in which a big event starts to quickly 
propagates through the whole cell in an almost periodic way; (c) nearly periodic regimes of 
a soliton-like character and (d) a regimes involving events with many size and time scales. 
In the latter case, they found exponential statistical distributions in amplitude, duration 
and separation time between events. 

With respect to theoretical and numerical studies, one of the systems used to model the 
dynamics of earthquakes has been the model introduced by Burridge and Knopoff (BK) in 
1967, and investigated extensively in recent publications 1^,0]. Although most of the studies 
on the Burridge-Knopoff model have been concentrated on the parameter region where a 
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partial power-law (partial here means of limited size) of the distributions of event sizes is 
found, this regime in fact has not been seen yet in experimental studies of a homogeneous 
system. The three first regimes (a), (b) and (c) found by Rubio and Galeano have been 
observed in numerical studies of the BK model. No reports so far exist of exponential 
distributions in the BK model. 

In this paper we show that exponential distributions are also seen in the BK model. This 
kind of distribution is observed when the two kinds of springs in the model have the same, 
or approximately the same, value for the elastic constants. They are present in quantities 
such as the time duration, number of blocks displaced, and total displacement (moment) of 
the blocks involved in an event. 

The paper is organized as follows: in the next section we describe the model; in the third 
section we study the statistical distributions and the different regimes that appear when the 
spring constants are varied; the last section is dedicated to the conclusions. 

II. THE BURRIDGE-KNOPOFF MODEL 

The homogeneous version of the Burridge and Knopoff model is shown in Fig. 1 of 
Ref. [0. It consists of a one dimensional array of blocks, each of mass m, coupled by 
springs of constant kc to one another, and by a spring of constant kp to a rigid pulling bar 
that moves at constant velocity V. In equilibrium, when all the springs are unstretched, 
adjacent blocks are separated by a distance a. The blocks rest upon a stationary surface, 
which provides a frictional force that impedes the motion of the blocks. In the version 
considered in Ref. the friction is a decreasing function of the velocity, the same for all 
blocks. The equation of motion for the jth mass when it is moving is 

mXj = kc{Xj+i - 2Xj + Xj^i) - kp{Xj - Vt) - F{Xj/Vf), Xj ^ 0, (1) 

where Xj denotes the displacement of the block measured with respect to the position where 
the sum of the elastic forces on it is zero. The last term in Eq. (|I|) represents the nonlinear 
and velocity dependent frictional force, which is given by 
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F(X,/V» 



Fo/(l + X,/V7), if X,>0 
-cx), if Xj < 

where Vj is the characteristic velocity for the friction. Here, we do not allow backward 
motions and the friction force gets arbitrarily high to avoid it. This is just for computational 
convenience, and this does not change the main results presented here. 

As shown in Ref. [0], this model has five independent velocity scales. Two appear ex- 
plicitly in Eq. (|l]): (i) the pulling velocity V and (ii) the characteristic friction Vj. Two 



more can be defined in terms of the spring constants, the mass and Fq: (iii) Vq = Fq/ Jmk 



and (iv) V/ = {F^/kp)^ kc/m. Here, Vq and Vi correspond to the maximum velocities of a 
single block held by a spring of constant kp and fee, respectively, when it has displaced by 
the characteristic distance Fo/kp in the absence of friction. The fifth velocity is the sound 



velocity Vg = a^kc/m, which depends on the equilibrium spacing of the blocks a and does 
not appear explicitly in Eq. (|l]). Depending on the relative sizes of these velocities, we would 
expect different behavior for the system. 

Following Ref. 0, we introduce different dimensionless variables 

U, = I^X,, r = {m/kpy/\ (2) 
-To 

so that Eq. can be written in dimensionless form 

t), = - 2f/, + - U, + ur- $(f7,/z//), (3) 



where 



1/(1 + U,/uf), if f7,>0 
-oo, if ilj < 



Now, z// = \JkJkp, V = V/Vq, and Uf = Vj/Vq. Dots here represent derivatives with respect 



to the scaled time r. The sound velocity becomes = a-'ykc/kp = a'ui, where a' = akp/kc- 
Without losing generality we can take the dimensionless length a' = 1. Thus, the sound 
velocity becomes identical to ui. The five velocities V, Vf, V/, Vg and Vq have been 



transformed respectively into u, Uf, ui, Vg = z//, and z/q = 1. Consequently, the model has 
three relevant parameters, z/, Uf and ui = Ug, which completely determine the behavior of 
the system. In the case of open boundary condition, which is the one we consider here, we 
have Uq = Ui and Un+i = Un- 

For fixed vi and v (large vi and small z/), a transition was reported in Ref. |^ when 
two velocities cross, namely, when Vf = vq = 1. For vj = oo one finds that the motion of 
the system is continuous. No block ever stops. As vj decreases, one sees small regions of 
stationary blocks. When Vf becomes less than 1 these stationary (event-free) regions begin 
to percolate across the entire system. The motion of the fault in this small vj region now 
occurs in abrupt large events. 

Here we report another transition, which occurs when the the velocity vi is varied for 
fixed values of v and vj. We find that for vi less than a given value the event sizes and event 
durations obey exponential distributions. Another value is observed for z/;, above which 
the distributions present partial power-laws, namely, events that have size smaller than a 
critical value have power-law distributions, and events that are larger than this value obey 
different statistics. In the intermediate regime for z/;, the distributions are neither power- 
laws nor exponentials. As vi gets close to zero, the exponential distributions disappear, since 
in that limit the blocks are disconnected, and only one-block events are observed. Evidence 
of a different behavior when vi is varied was found in Ref. [§. However, no statistical 
distributions were calculated in that paper. 

In the numerical simulations of this letter we have generally started the system with the 
blocks at rest and with the sum of the elastic forces in each block equal to zero. Other 
initial conditions were considered and the results did not change. We have also considered 
periodic boundary conditions, and the results were the same (except for the few events 
involving the boundaries of the chain). Before we start to compile statistics we let the 
system evolve until it reaches a statistical stationary state. Here we fix the pulling velocity 
and the characteristic velocity to the following values: v = 0.01 and uj = 1/6. It is beyond 
the scope of this paper to do a detailed study on the effects of variations in the pulling 
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velocity and in the characteristic speed. Studies on event distributions for varying z// and 
V with fixed vi can be found in [6]. We expect that the kind of transition we report here 
will be observed for other values of Vf and z/, and the value of vi where the transition occurs 
probably varies with the characteristic and pulling velocities. 

III. NUMERICAL RESULTS 

To illustrate the different behaviors of the system when vi is varied, we show in Fig. [l| 
projections of the block velocities Vj onto the j — t plane. A black dot in the figure means 
that the j-block is moving at time r and a white one, that it is at rest. We show three 
different cases with ui = 1,3,10. It is clear that a transition occurs as z/; is changed. For 
small i>i, we see clusters of moving blocks, which involve only a small number of masses. 
As ui is increased, one finds clusters involving a small number of blocks, as well as clusters 
involving a large number of masses. These large ruptures appear V shaped and travel with 
the sound speed z/;. 

Now we investigate the distributions of moment and time duration for different values 
of i>i. The moment of an earthquake is a measure of its size, and in dimensionless units it 
is defined as m = J2j ^Uj, where 6Uj is the displacement of the j-th block and the sum is 
over the blocks displaced during the event. The distribution r(m) is the number of events 
with magnitude m divided by the number of blocks in the chain and by the total time of 
dynamical evolution. For z/^ < 2 we find that r(m) is governed by an exponential function, 
namely, r(m) ~ exp(m/mc), where rric is a characteristic moment. The result of r(m) with 
ui = 1 is shown in Fig. ||( solid line. We have also found that the characteristic 

moment decreases as ui becomes smaller, since when 1^1 = the blocks become disconnected 
and only events involving one block are observed. For intermediate values of ui (2 ^ z/; ^ 4) 
ones sees a regime for r(m) which is neither power-law nor exponential behavior. This is 
shown in Fig. |^(b) where the solid line represents r(m) for Ui = 3. For large Ui (yi ^ 4) 
one sees the appearance of a partial power-law distribution. This was the parameter region 
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investigated by Carlson and Langer 0. Beyond the power-law cutoff one sees a bump in the 
distributions, which increases in height as ui becomes larger. We have found that the events 
that do not belong to the scaling region become more and more frequent as ui is increased. 
Thus, in the limit of very large ui the power-law tends to disappear. In Fig. ^(c) the solid 
line represents r(m) for Ui = 10. 

For the distributions of time durations 6t of the events we find qualitatively the same 
kind of behavior as found for the distributions of moments. For z// ^ 2 we get an exponential 
function, r((5r) ~ exp{ST / 6tc) . This is shown in Fig. ^(a) for z/^ = 1 as a dashed line. Again, 
we have found that the characteristic time duration 6tc decreases as ui gets smaller. For 
z// > 4 one sees partial power-law distributions, as shown in Fig. ||(c) for z/^ = 10 as a dashed 
line. For intermediate values of z// (2 ^ z/; ^ 4) the distributions are neither exponential nor 
power-law ones, as we show in Fig. ^(b) for z/; = 3 (dashed line). 

We have also studied the statistical distributions r(n) of an event involving n blocks. 
Again, the same kind of behavior described above is seen, namely, exponential distributions 
for ui ^ 2, partial power-laws for ui ^ 4, and a regime without any scaling features for 
intermediate values of z/;. 

IV. CONCLUSIONS 

We have found exponential distributions for the sizes and time durations of the events in 
the Burridge-Knopoff model for earthquakes when the two kinds of springs have the same 
(or approximately the same) elastic constants. In this situation the velocity ui and the 
sound velocity are equal (or approximately equal) to one. For intermediate values of ui, the 
distributions are neither exponential nor power-law. For large ui, partial power-laws (namely, 
a scaling region of limited size) are observed. Since in the experimental system studied by 
Rubio and Galeano the elastic gel is homogeneous, this implies that the exponential 
distributions they found are consistent with our results for the case in which the springs of 
the BK model have identical elastic constants. 
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FIGURES 

FIG. 1. Projections of the block velocities Uj onto the j — r plane for a 200 block system with 
u = 0.01, Uf = 1/6, and (a) = 1, (b) ui = 3, and (c) = 10. 

FIG. 2. Statistical distribution r{m) of the moment of the events (solid lines) and the distri- 
bution r((5r) of the time durations of the events (dashed lines) for (a) ui = 1, (h) vi = 3 and (c) 
ui = 10. In all cases, = 300, u = 0.01 and Vf = 1/6. 
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FIG. 1(a) 
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FIG. 1(b) 
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FIG. 1(c) 
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